Frontiers in Physics

The COVID-19 pandemic triggered an unprecedented volume of real-time discourse on social media platforms, with Twitter serving as a global forum for public reactions, fears, and evolving narratives. Traditional sentiment analysis approaches treat tweets as independent, static samples, failing to capture the temporal evolution and geographic heterogeneity of public opinion. This paper presents a comprehensive spatio-temporal framework that integrates fine-grained sentiment classification using COVID-Twitter-BERT with dynamic topic modeling via BERTopic to automatically discover and track evolving narratives. Using a corpus of 2.4 million geolocated tweets collected between January 2020 and June 2022, our analysis reveals distinct pandemic phases: early fear-driven narratives about mask shortages (Q1 2020), vaccine optimism followed by polarization (2021), and pandemic fatigue (2022). Regional comparisons show significant differences, with US discourse dominated by freedom-versus-mandate debates while European discussions emphasized collective solidarity. Our framework achieved 76% F1-score in sentiment classification and successfully identified 50 distinct narratives with high coherence scores. This work provides a powerful methodology for real-time epidemiological narrative surveillance and crisis communication monitoring.

A fundamental biological process of transcription occurs in the cell nucleus, which is a complex medium that also contains multiple heterogeneous structures known as biomolecular condensates. Interestingly, some of these condensates contain mRNA molecules in addition to proteins, suggesting an important cellular role in transcription that is not yet well understood. In this work, we develop a minimal theoretical framework for quantitative investigation of the role of reversible mRNA condensation in transcription. Our discrete-state stochastic approach accounts for the most relevant processes, allowing us to explicitly evaluate the properties of the system and clarify the effects of condensation. Analytical calculations supported by computer simulations suggest that reversible mRNA condensation influences the transcription processes by maintaining a constant level of free mRNA in the nucleoplasm while lowering the degree of stochastic noise and increasing the robustness against external perturbations. Physicochemical arguments are presented to explain these observations. The proposed theoretical framework elucidates important microscopic aspects of transcription, providing a convenient quantitative tool for investigating complex biological phenomena.

In the face of an epidemic where a population behaviour both influences disease transmission and reacts to it, social processes can generate norms to support socially beneficial behaviour. Most mathematical models of coupled behaviour disease dynamics treat norms as pre-existing rather than explaining how they are maintained. Here, we investigate whether altruistic punishment can sustain a social distancing norm when individuals may defect, cooperate without punishing, or cooperate while paying a cost to punish defectors. We couple a transmission model to an imitation model for these three strategies. Disease prevalence affects behavioural payoffs, while the behavioural composition modifies transmission. We also compare this coupled system with a control where behavioural decisions respond to a fixed prevalence. We find a wide parameter regime corresponding to the establishment of an injunctive social norm in support of social distancing, where the punisher strategy is widespread. Persistence may occur through stable states where punishers or dominant. Disease behaviour feedback can also create oscillations (where the three strategies succeed one another in response to epidemic waves) or tipping points (sharp transitions between all-defector and cooperative states). These effects do not occur in the uncoupled model, although there are still broad parameter regimes where a social norm persists. Our findings show that costly peer punishment can support persistence of social norms that mitigate disease transmission. More broadly, endogenous epidemic feedback can qualitatively change the conditions under which cooperation and punishment are sustained, producing tipping points and long-term behavioural epidemiological cycles that fixed-payoff models cannot capture.

Biomolecular condensates are widespread cellular self-assembled structures with essential functions. There are suggestions of condensates formed by different proteins being near criticality. However, systematic investigation of the criticality of condensates is absent, and critical exponents defining their universality class have not been found. Here, using long-time simulations, we show that condensates exhibit typical critical phenomena, including scale-free spatiotemporal correlations, critical slowing down, divergence of correlation length and dynamic scaling. From these scaling behaviors, a set of critical exponents is determined. Based on dynamic critical exponent, diverse condensates can be divided into two distinct universality classes, arising from differences in their molecular components and interaction types.

Mechanical properties of biological tissues, driven by passive and active forces, play a vital role in several processes ranging from development to cancer metastasis. However, the dynamical responses of cells in tissues, subject to mechanical deformations such as shear and the associated rheological properties, are not well characterized. Here, we use three-dimensional agent-based models for normal and cancer tissues to investigate their responses to simple shear as a function of cell stiffness and stochastic active forces. In the normal epithelium, with uniform strength of active force, the yield stress as a function of shear rate follows the Herschel-Bulkley form over a range of cell volume fraction. Strikingly, the shear rate dependence and the elasticity-dependent changes in the yield stress fall on master curves upon suitable scaling. To model cancer-like behavior, a certain fraction (Np) of cells was chosen to have enhanced activity and decreased stiffness. As Np increases, the extent of collective cell movement decreases, transitioning from affine (collective) to non-affine (individualistic) movement, a finding that is in accord with imaging experiments. Simulations of a model of a stiff solid tumor, with radius Rs embedded in normal tissue, show that as Rs increases, the yield stress increases. Interestingly, the cells migrate collectively as Rs increases. A Gaussian Mixture Model (GMM) and a mean field theory quantitatively account for the simulation as well as experimental results on cancerous, non-cancerous, and a mixture of these two types. The combined theoretical and experimental study establishes that heterogeneity in stiffness and activity determines non-affine movements in normal and cancer tissues.

Cardiac blood flow is a regulator of several important developmental and remodelling processes in the heart, including through fluid shear forces sensed by the endothelial cells lining the heart. However, optically mapping these flow fields in the complex 3D geometry of the heart is challenging even in transparent animal models such as the zebrafish. One of the main challenges is the difficulty in measuring the out-of-plane (axial) velocity component, preventing accurate mapping of the complete 3-component-3-dimension (3C-3D) blood flow velocity field; image-based techniques such as microscopic particle image velocimetry ({micro}PIV) traditionally only provide the in-plane flow components. Here we present a computational approach to achieve full time-varying 3C-3D blood flow vector mapping using a standard selective plane illumination microscope (SPIM), based on robust cardiac phase assignment, precise measurement-driven registration of sequentially acquired z-stacks, and PIV data fusion from multiple sample orientations. Our approach holds the key to understanding the complex dynamic flow fields within the developing heart, and their role in shaping cardiac development.

AO_SCPLOWBSTRACTC_SCPLOWCellular automata and graph reaction-diffusion systems encode local spatial interactions in different mathematical forms. We develop a cochain-operator calculus for these two settings. Over a finite field Fq, every local rule on a finite neighborhood has a unique reduced polynomial representative. On an oriented line, the coboundary and endpoint maps recover the left and right shifts. Our main theorem shows that these operators, together with linear operations, constant cochains, and the degree-zero cup product, generate every finite-radius polynomial cellular automaton. Explicit formulas for Rules 30, 110, and 22 show how reflection-invariant linear coupling, directed transport, and nonlinear neighbor interactions enter the calculus. On a general graph, d*d is the unweighted combinatorial Laplacian and enters a graph reaction- diffusion recurrence. Over [R], the term - Dd*d with D [≥] 0 admits the usual diffusion interpretation; over Fq, the corresponding expression defines modular coupling without an intrinsic order. In the morphogenetic examples, we therefore distinguish pattern-generating dynamics from finite-state observation and use the Betti numbers of active induced subcomplexes to summarize observed patterns. This yields a common algebraic representation without identifying real-valued diffusion with finite-field dynamics.

We develop a theoretical proposal linking vacuum stability and brain dynamics through superconductivity-inspired coherence, symmetry reduction, and the thermodynamic stabilization of low-entropy regimes. We take an unbroken SU(3) structure as a candidate stable residue of the low-temperature vacuum. At the neural level, we formulate a coarse-grained analog in which a two-fluid model with dissipative and coherence-supporting components describes brain dynamics. Specifically, the coherence-supporting component is proposed as a possible basis for the efficient binding and integration required to sustain a stable, unified conscious state. The proposal offers a common geometric language for relating physics and neuroscience with falsifiable signatures in coherence and state-dependent transitions. The main technical contribution is a computational algebraic model of conscious-state dynamics, where neural data are mapped to reconstructed state trajectories. Effective generators are inferred from those trajectories, and the two-fluid split is tested as a Cartan-root decomposition of su(3), with a rank-two commuting sector for coherence-preserving balance and six root directions for state transitions. This structure can be tested on neural data and contrasted with alternative dynamical models.

The coexistence of competing pathogen strains is shaped by cross-immunity, the cross-protection that infection with one strain confers against another. Although cross-immunity is often asymmetric between strains, this asymmetry is often neglected in the literature on multi-strain coexistence. The effect on coexistence-exclusion outcomes of waning immunity\textemdash which is particularly relevant for antigenically evolving pathogens\textemdash is also poorly understood. To understand how these factors affect strain coexistence, here we analyze a status-based two-strain SIRS model with asymmetric cross-immunity and strain-specific rates for transmission, recovery, and waning of immunity. We derive explicit invasion thresholds that also determine the feasibility and local stability of a unique coexistence equilibrium. Thus, these thresholds allow us to characterize the region of stable strain coexistence, as a function of the cross-immunities and rates of waning immunity. We also obtain closed-form expressions for the strain prevalences at the coexistence equilibrium, showing that the total prevalence may vary non-monotonically as the basic reproduction number of one strain increases. Finally, we show that a transient reduction in transmission can move a coexisting strain pair across an invasion boundary, driving the weaker strain extinct. Applying this result to influenza B, our analysis offers a parsimonious explanation for the disappearance of the Yamagata lineage during the COVID-19 pandemic.

Interleukin-13 (IL-13) is an immunomodulatory cell signaling cytokine that has been implicated in neurodegenerative disease and chronic inflammation. IL-13 binds to its low and high affinity receptors, IL-13 receptor 1 (IL-13R1) and IL-13 receptor 2 (IL-13R2), respectively, with residence times that vary accordingly. As the binding kinetics of the cytokine-receptor complexes influence cellular responses, we employed the molecular dynamics (MD) simulation-based{tau} -random acceleration molecular dynamics method ({tau}RAMD) to compute relative residence times for wild-type (WT) IL-13 and 19 IL-13 mutants to the two receptors. Comparison with experimental kinetic data shows that the{tau} RAMD computations capture the trends in residence times. Analysis of simulated dissociation trajectories of the cytokine-receptor complexes reveals two distinct dissociation pathways of IL-13 from each of the receptors. This study thus pinpoints key determinants of the interaction of IL-13 with its receptors which could be targeted for therapeutic design. Statement of SignificanceCytokines are regulatory proteins that bind to cell surface receptors and thereby send signals to the cellular interior. Interleukin-13 (IL-13) is a cytokine that has a low and a high affinity receptor. It has important physiological roles, and its deregulation is involved in diseases such as atopic dermatitis and asthma. We employed a molecular dynamics simulation-based method to compute the effects of changes in the sequence of IL-13 on the lifetimes of complexes of IL-13 and its receptors. Comparison with experiments supports the validity of the computational approach and analysis of the simulations reveals two distinct ways in which IL-13 dissociates from each receptor. These results thus provide a map for targeting IL-13 - receptor interactions for the design of therapeutics.

Biomolecular condensates are dynamic, membrane-free compartments that continuously exchange molecules with their surroundings. The dwell time, defined as the time a molecule remains inside a condensate between entry and exit, determines how extensively the molecule can explore the dense phase and encounter potential binding partners or reaction sites, thereby modulating condensate function. Motivated by our single-molecule measurements of RNA dwell times, we developed an analytical theory to understand dwell-time distributions in biomolecular condensates. Our theory predicts that the dwell-time distributions generally exhibit an early-time power-law regime followed by a late-time exponential tail. The form of the distribution encodes the rate-limiting mechanism of molecular escape: dense-phase diffusion-limited transport feature a -1.5 power law with an exponential tail set by a diffusion timescale, whereas interfacial barrier-crossing-limited transport feature a -0.5 power law with a decay governed by a barrier-crossing timescale. These distinct signatures provide a direct readout of the physical processes that control molecular retention in condensates, with implications for both natural and synthetic condensates.

Wound repair is an important mechanism to preserve tissue integrity in organisms after injury. However, why different tissues exhibit different mechanisms to repair wounds is a long-standing question that remains unanswered. In this work, we theoretically explore the role of the purse string, an actomyosin contractile cable used by tissues to close small wounds. Does the tissue 3D geometry influence the efficiency of the purse string in driving wound closure? Using a 3D biophysical model, we study in silico tissues with the same cell volumes but different aspect ratios, ranging from squamous to thick and tall tissues. The model predicts that taller cells are easily deformed by the purse string. In contrast, very squamous cells require a very strong purse string that might demand additional cellular mechanisms to close the gap. These findings establish a theoretical framework to predict the optimal biophysical mechanisms of wound healing in different tissues. Graphical abstractCells of different aspect ratios can be observed in a range of organisms with different function and mechanics. The wound healing efficiency of the purse string increases with the cell aspect ratio in our theoretical exploration. O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=135 SRC="FIGDIR/small/743165v1_ufig1.gif" ALT="Figure 1"> View larger version (23K): org.highwire.dtl.DTLVardef@d44ab0org.highwire.dtl.DTLVardef@1737cbaorg.highwire.dtl.DTLVardef@101b5d4org.highwire.dtl.DTLVardef@1487f26_HPS_FORMAT_FIGEXP M_FIG C_FIG

In a population model for an infectious disease, we consider the early stochastic dynamics of an emergent 'mutant' strain, appearing and spreading during an epidemic of another 'wildtype' strain. The mutant may not reach establishment in the host population. The time at which the mutant first appears determines its probability of establishment. We calculate this establishment probability with two methods. The first method assumes a classical branching process, with a constant transmission rate. The second method reflects the changing size of the pool of susceptible hosts, due to the dynamics of the wildtype. We find that susceptible depletion can substantially impact the establishment probability. We explore the consequences of this stochastic establishment on the "escape pressure" acting on a pathogen to produce immune escape variants. We find that the overall escape pressure rate depends strongly on the appearance time of the mutant, especially if the establishment probability is itself shaped by the continued spread of the wildtype. In most scenarios, the escape pressure rate (and thus, the risk of new escape variants) peaks slightly earlier than the prevalence of the wildtype strain. Integrating the escape pressure over time, we obtain the cumulative escape pressure generated by the wildtype epidemic. The relationship between the escape pressure and the vaccination coverage depends on the cross-immunity, due to susceptible depletion. For example, with intermediate cross-immunity, the risk of immune escape may be lowest at intermediate vaccination coverages. Thus, these results raise important considerations for vaccination strategies in response to novel outbreaks.

In the active loop extrusion model, the cohesin protein complex creates chromatin loops in eukaryotic cells. Extrusion maintains topologically associated domains (TADs), which are contiguous segments of chromatin that preferentially colocalize in space and are typically bounded by CTCF proteins that pause cohesin translocation. Here, we model active loop extrusion with hybrid molecular dynamics - Monte Carlo simulations in entangled flexible linear polymer melts. Intra-chain contact probabilities of polymers with active loop extrusion are enhanced compared to their equilibrium, passive counterparts. Extrusion causes the size of chain segments to be much smaller than in passive melts. While the overlap parameter in passive melts without extrusion monotonically increases with segment length, it is nonmonotonic in active melts and on the order of unity within the parameters of this study. Active loop extrusion suppresses contacts between TADs in favor of intra-TAD contacts. Reduction of overlaps between chain segments dilutes entanglements in active melts. Depending on parameters, active extrusion without TADs may induce more compact conformations than with TADs, due in part to fractal loopy globule-like dynamics. This work suggests that active loop extrusion reduces overlaps between TADs, contributing to effective gene regulation by cis-regulatory elements.

Differences in the density of heterochromatic and euchromatic regions are often assumed to play an important role in determining the transcriptional state of chromatin by influencing its accessibility to transcription factors. Yet, experiments show that even the most densely packed chromatin domains are readily accessible to transcription factors. Thus, the molecular mechanisms underlying differences in the transcriptional states of hetero- and euchromatin remain poorly understood. In this study, using electrically charged mEGFP probes, we demonstrated that heterochromatin and euchromatin differ not only in density but also in the magnitude of the local electrostatic potential. Furthermore, the nucleoplasmic distribution of electrically charged proteins, such as histone-chaperone complexes, was found to correlate with the local electrostatic potential. Estimates based on experimental data have also shown that the difference in electrostatic potentials between heterochromatic and euchromatic regions could lead to unequal nucleosome stability in them, which was successfully confirmed experimentally. Subsequent theoretical calculations showed that this could shift the balance in the DNA-binding competition between transcription factors and nucleosomes, thereby explaining experimental observations that heterochromatin is less transcriptionally active than euchromatin. Overall, our study suggests the existence of a nuclear electrostatic potential-mediated pathway that may be involved in the regulation of gene transcription.

To disentangle complex fiber configurations that remain challenging for diffusion MRI tractography, insights might be gained from microscopy tractography. Indeed, by precisely following small white matter (WM) fascicles invisible at the resolution of diffusion MRI, microscopy tractography can help explain how fiber populations are organized at the finest scales. Serial optical coherence tomography (S-OCT) is an imaging modality relying on the intrinsic contrast of a sample. When applied to brain tissues, the S-OCT contrast is primarily driven by the myelin reflectivity. Due to its high resolution, on the order of microns, and its 3D nature, S-OCT offers promise for studying WM connections at the microscale. However, while other microscopy imaging modalities have been shown to enable tractography, whether the reflectivity contrast from S-OCT supports the reconstruction of long-range WM fascicles at the microscale remains unknown. Furthermore, there is a gap in the literature regarding how an ideal microscopy tractography algorithm should behave with respect to the choice of tractography algorithm, tracking maps definition and microscale orientation distribution functions (ODF) estimation. In this work, we describe a tailored approach to reconstruct WM fascicles at the microscale from S-OCT acquisitions. We improve microscale orientation distribution functions (ODF) estimation by implementing a sliding-window formulation allowing the estimation of ODF at S-OCT resolution, and use apodized Dirac delta functions for reducing unwanted interference. We validate our approach on a simulated microscopy-like FiberCup dataset, and show that using multiscale Frangi filters for estimating ODF outperforms structure tensor analysis. We also show that particle filtering tractography with anatomical constraints enables targetted, region-to-region tractography, and outperforms standard deterministic or probabilistic tracking approaches. We further demonstrate our method on a whole mouse brain S-OCT reconstruction at 10 m by reconstructing the thalamocortical white-matter projections. Overall, our results show that S-OCT tractography recovers fine white matter fascicles visible at the microscale, and that these connections are supported by viral tracing experiments from the Allen Mouse Brain Connectivity Atlas. Moreover, this work shows the first ODF estimation and fully-3D probabilistic particle filtering tractography of the mouse brain from S-OCT reconstructions at 10 m isotropic resolution.

Living cells are highly dynamic and densely crowded environments in which organelles such as vesicles undergo continuous motion that is essential for cellular processes. Therefore, accurate tracking of individual organelles is crucial for understanding intercellular dynamics and functions. However, precise tracking of individual organelles in living cells remains challenging due to high organelle densities, frequent particle overlap, and the coexistence of stationary and motile organelles. In particular, stationary organelles can obscure the trajectories of moving organelles, leading to tracking errors and fragmented tracks. To overcome these challenges, we developed Multiple Particle Tracking via Velocity Filtering (MPT-vVF), an unbiased, semi-automated tracking framework that incorporates a mathematically derived velocity-filtering algorithm to selectively identify and track moving organelles with high accuracy in crowded intracellular environments. MPT-vVF integrates denoising, background subtraction, and a velocity-matching detection step that discriminates true particle motion from noise based on spatiotemporal continuity, followed by robust trajectory linking. We demonstrate that MPT-vVF can accurately resolve nanometer-scale displacements of immobilized beads, highlighting its high tracking precision. We also validate the robustness of MPT-vVF by quantifying the transport of brain-derived neurotrophic factor (BDNF)-mRFP-containing vesicles in living hippocampal neurons. Furthermore, MPT-vVF reveals that exposure to 50-nm nanoplastics impairs vesicular transport, reducing both travel length and speed of BDNF-containing vesicles in living neurons. These findings establish MPT-vVF as a powerful method for quantitative analysis of intracellular organelles in crowded living cells and suggest its broad application to biophysics, cell biology, and soft matter research.

Biomolecular condensates are dynamic assemblies whose functions depend on continuous exchange of molecular components with the surrounding environment. While scaffold molecules drive phase separation and condensate architecture, many functional components are clients that are recruited through interactions with the scaffold-rich environment. Despite their prevalence, how client-scaffold interactions shape client exchange dynamics remains poorly understood. Here, we develop a reaction-diffusion model for client exchange in scaffold-driven condensates, in which clients switch between a scaffold-bound state and an unbound state. Bound clients exchange through scaffold-mediated transport, whereas unbound clients diffuse through the pore space of the condensate. Using the fluorescence recovery of fully photobleached condensates as a measure of client exchange, we compare transport through these two pathways with bound-unbound conversion and identify three limiting regimes. In the slow-conversion regime, bound and unbound clients recover through distinct scaffold- and pore-mediated pathways. In the intermediate-conversion regime, recovery of bound clients becomes limited by client unbinding. In the fast-conversion regime, local equilibrium between bound and unbound clients produces an effective single-state recovery. We further propose a unifying description that connects these regimes and quantitatively captures the apparent recovery timescales extracted from numerical simulations across condensate sizes. Our results provide a framework for interpreting component-specific exchange dynamics, and highlight client size, client-scaffold binding, and condensate porosity as key regulators of client turnover in multicomponent condensates.

Protein folding is the process by which a polypeptide chain organizes into its three-dimensional structure through a balance of stabilizing and destabilizing interactions encoded by the sequence. A central question in protein biophysics is how thermodynamic factors guide a polypeptide toward its native folded state despite the rugged energy landscape and the competing influence of nonnative interactions. In many biomolecular processes, cooperativity provides a mechanism by which multiple weak interactions act collectively to generate a robust response. In the context of protein folding, such cooperative effects may arise when the formation of one native contact enhances the stability or likelihood of nearby native contacts, thereby promoting collective organization toward the folded state. At the same time, folding is opposed by the much larger number of non-native interactions, whose heterogeneity can introduce frustration and destabilize folding even when the average native bias favors the folded phase. The interplay of these competing effects in determining foldability remains unclear in statistical-mechanical models. Here, we address this problem using a one-dimensional spin-glass model of protein folding with explicit shared-residue cooperative interactions encoded through wedge-based motifs. We show that modest cooperative bias can stabilize folding even where the noncooperative system remains unfolded, whereas non-native energetic fluctuation suppresses folding and shifts the transition to higher cooperative strengths. We further find that partial cooperative coverage is sufficient to lower the folding threshold. Therefore, the model provides a mean-field framework for incorporating cooperative interaction strength into the native one-dimensional model of protein folding and for describing how local cooperativity reshapes the folding transition.

While geometric constraints shape tissue development, quantifying the resulting growth dynamics remains a central challenge in tissue engineering. Conventional methods often struggle to capture multi-scale kinetics without complex labeling or difficult single-cell tracking. Here, we analyze geometrically controlled growth of microtissues derived from human dermal fibroblasts using time-resolved, label-free brightfield microscopy, combined with optical flow and semi-automated deep learning mitosis detection. By extracting multi-scale flow fields and integrating them with tissue segmentation, we quantify directional tissue dynamics, separating flow into parallel and normal components relative to the local tissue contour. Applying this framework, we contrast the quiescent tissue interior with the advancing growth front where localized dynamics and cell proliferation drive expansion. Our results demonstrate that, compared to the bulk, the growth front exhibits higher fluctuations parallel to the tissue contour, positive mean normal flow, and significantly increased mitotic activity. Furthermore, evaluating flow divergence around mitotic events reveals distinct spatial behaviors: with the onset of mitosis, a contraction and subsequent expansion occurs in the vicinity of the dividing cells. Beyond the immediate cellular neighborhood, the broader regional dynamics remain consistent before and after mitosis onset, with net tissue expansion in proximity to the growth front and contraction within the tissue interior. By extracting continuous kinetic data from easily accessible, label-free brightfield imaging, this approach serves as a non-invasive, complementary tool for evaluating in vitro tissue morphogenesis and growth dynamics. This analytical framework can be expanded to study locally resolved tissue morphogenesis and growth kinetics in other microsystems, ranging from embryos to organoids. Statement of SignificanceUnderstanding how localized cellular forces drive tissue growth is critical for mechanobiology. However, mapping these dynamics traditionally requires complex, invasive fluorescent labeling. We present an accessible, label-free computational framework combining optical flow and deep learning-based mitosis detection to quantify continuous tissue kinematics directly from standard brightfield microscopy. Applying this to 3D microtissues, we reveal a distinct spatial coupling between cell division, local mechanical fluctuations, and directed tissue expansion at the active growth front. This non-invasive approach bridges the gap between single-cell mechanics and macroscopic morphogenesis, offering a versatile tool to monitor complex in vitro model systems-like organoids and bioengineered tissues-without disrupting their native state.